Non-standard pairing in asymmetric trapped Fermi gases 
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We study an ultracold trapped Fermi gas of atoms in two hyperfine states with unequal popula- 
tions. In this situation the usual BCS pairing is suppressed and non-standard pairing mechanisms 
become important. These are treated by solving the Bogoliubov-de Gennes equations, which at the 
same time correctly take into account the finite size of the trapped system. We find results which 
can be viewed as generalization of the LOFF phase to finite systems. 
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In the last years much progress has been made in im- 
proving the techniques used to trap and cool dilute gases 
of bosonic and fermionic atoms 0,0- One of the interest- 
ing aspects of the properties of ultra-cold gases is that the 
interatomic interaction can be modified, both in its in- 
tensity and in its sign, by changing the applied magnetic 
field around a Feshbach resonance. Due to the very low 
densities and temperatures in these systems the details of 
the interatomic interaction are unimportant and the in- 
teraction can be characterized by one single parameter, 
the s-wave scattering length a. In this article we con- 
sider Fermi gases trapped and cooled in two hyperfine 
states with an attractive interaction, i.e., a < 0. We will 
concentrate on the weakly interacting case (fcp|a| <C 1, 
where kp denotes the Fermi momentum). In this region, 
BCS superfluidity with formation of Cooper pairs is ex- 
pected below a certain critical temperature. So far, some 
experimental signals have been found which would indi- 
cate the existence of superfluidity in this region [2j, but 
a clear evidence is still missing. 

Besides the interaction, also the population of the 
two hyperfine states can experimentally be controlled. 
Usually mixtures with equal populations are created, 
but controlled asymmetric mixtures have also been ob- 
tained Q. Unequal populations of the two hyperfine 
states can lead to very interesting phenomena. For in- 
stance, the BCS pairing mechanism is supposed to be- 
come suppressed Q since the two Fermi momenta as- 
sociated with the two spin polarizations become differ- 
ent: The formation of zero- momentum Cooper pairs built 
of two atoms at their respective Fermi surface becomes 
difficult. Instead, other more exotic pairing phenomena 
have been suggested for the case of unequal populations, 
like the Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) phase 
0, the Sarma (interior gap) phase 13 Q, or a phase 
with deformed Fermi surfaces (DFS) J8j . Many of these 
non-standard pairing mechanisms have already been dis- 
cussed in other domains of physics where asymmetric 
two-component fermion systems can be found: Super- 
conductors in a magnetic field 0, IT(i| . n eutron-proton 
pairing in asymmetric nuclear matter color super- 
conductivity in high density QCD |12| . The experimen- 
tal observation of the LOFF phase in the case of super- 



conductivity is a subject of debate. It seems that quite 
recently an oscillating superconducting order parameter 
has been observed in a ferromagnetic thin film |l3[ . 

Usually 0, IE IE @ these non-standard pairing types 
in ultracold Fermi gases are usually discussed for the case 
of infinite and homogeneous systems, or for trapped sys- 
tems in local-density approximation (LDA) which again 
amounts to treating the system locally as homogeneous. 
However, as we are going to show, in certain cases the 
LDA can become a very poor approximation and we 
therefore want to analyze this problem within a fully 
microscopic mean field Bogoliubov-de Gennes (BdG) ap- 
proach |14| taking into account the inhomogeneity and 
finite size of the trapped system. Recently the solution 
of the BdG equations has been considered in Ref. [l5j . 
where the authors discuss also possibilities for the exper- 
imental detection of a spatially modulated pairing gap. 

In the present article we study two systems: a small 
one where shell effects still play a role and a relatively 
large one, where one could expect the LDA to be a rea- 
sonable approximation. As we will show, this is not al- 
ways the case, although the LDA describes roughly the 
qualitative features. In addition, we examine the tem- 
perature dependence of the non-standard pairing phase, 
since this is important in connection with the experimen- 
tal realization of such a phase. 

Let us briefly recall the BdG formalism. We consider 
a system containing N = N + + N- atoms of mass m 
trapped by a spherical harmonic potential in two hyper- 
fine states denoted by + and — . The many-body Hamil- 
tonian of the system can be written as 

+ .gf 1 f (r)* 1 L(r)*_(r)'F + (r)) , (1) 



where H = -fr 2 V 2 /(2m) 



2 /2 denotes the har- 



monic oscillator single-particle Hamiltonian, ^(r) is the 
field operator which annihilates a particle at the position 
r in the spin state a (+ or — ), g = 4Trh 2 a/m is the 
coupling constant. In mean-field approximation, on can 
derive the following BdG equations corresponding to the 
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Hamiltonian QJ: 
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1<T (r)E ntr = W a u na {r) + A(r)« 7) _ CT (r) , 
^(r)^„ CT = -W_ CT v„_ CT (r) + A(r) V (r) , 



where W a = Hq + gp_ CT (r) — fi a , \i a ancl Pa being the 
chemical potential and the density, respectively. A(r) 
denotes the pairing field (gap) and E na , u va and Una arc 
the quasiparticle energy and wave functions, respectively, 
corresponding to the quantum numbers rf = n,l,m and 
spin a. In order to have different populations, the two 
chemical potentials fi + and fi- must be different. It is 
convenient to introduce the average chemical potential p 
and to write fi a = p + a8fi/2 where Sfi determines the 
asymmetry. Eqs. © reduce to the usual BdG equations 
in the symmetric case fi + = They arc solved numer- 
ically employing the same regularization method for the 
pairing field as described in Ref. |k| f° r the symmetric 
case. 

Eqs. are general enough to describe rather compli- 
cated types of non-standard pairing. In the case of usual 
BCS pairing, the dominant matrix elements of the pair- 
ing field are the diagonal ones, i.e., each state \n, I, m, +) 
is paired with its time-reversed counterpart |n, I, —m, —)■ 
However, the non-diagonal matrix elements of A are also 
included, which amounts to taking into account also the 
pairing of states \n,l,m,+) and n', I, —to, — } with differ- 
ent principal quantum numbers n ^ n. In our present 
calculation, we still keep the restriction that the Cooper 
pairs have total angular momentum zero. To release 
this constraint would mean to allow for a spontaneous 
breakdown of spherical symmetry, which would be nu- 
merically very heavy. The effect of strong non-diagonal 
matrix elements of A in fact corresponds closely to the 
LOFF phase in the case of a uniform system. There, the 
states are labeled by their momentum k. In the sim- 
plest version of the LOFF phase, the Cooper pairs have 
total momentum q, i.e., each state |k, +) is paired with 
— k + q, — }. The corresponding gap is oscillating with 
wave vector q and its matrix elements are therefore of 
the form A kk < = A<5 k _qj,< ( m contrast to the BCS phase, 
where A kk ' = A£ kk /) [l7j . 

In the discussion of our results all quantities will be 
given in harmonic oscillator units. We use the same cou- 
pling constant as in Ref. |l6j . i.e., g = —h^lho/m, where 
Iho = \/hJ{rrujj) denotes the harmonic oscillator length, 
and we consider two values for the average chemical po- 
tential, p = 22 Hui (small system with N w 4900) and 32 
hid (large system with N « 17000). 

In Fig. 1 we show the values of the pairing gap A(0) 
at the center of the trap for increasing aymmetry Sfi at 
T = 0. Let us first look at the lowest line, corresponding 
to the small system with p = 22 Hui. When both spin 
states are equally populated (Sfi — 0), we find A(0) « 2 
hw, i.e., we are no more in the intrashell-pairing regime, 
but shell effects are still important ^(|. If we increase Sfi, 
the two Fermi surfaces become more and more separated, 
i.e., if the state \n, I, m, +) lies close to the Fermi level for 




FIG. 1: Value of the pairing gap at the center of the trap 
(in units of hui) as a function of the asymmetry Sfi (in units 
of hu)). The lowest line corresponds to fi — 22ftu. The two 
upper lines correspond to fi = 32hu and show the BdG (solid 
line) and the LDA result (dashed line). 



spin +, the state \n, I, — m, — ) lies far from the Fermi level 
for spin — , making BCS pairing less and less favorable. 
As a consequence, at Sfi « 1.2 Hu>, corresponding to a 
particle number asymmetry a = (N+ — N-)/N ~ 0.06, 
the pairing disappears (shell closure effect). But then, 
near Sfi « 2 hw (a w 0.07), the states \n, I, m, +) near the 
Fermi level for spin + approach the states \n', I, —m, — ) 
near the Fermi level for spin — if n' = n — 1. Therefore, 
pairing becomes again possible, but now the Cooper pairs 
are built of two wavefunctions with different numbers of 
nodes, leading to a gap A(r) which as a function of r has 
exactly one node. 

Let us now turn to the investigation of the larger sys- 
tem, fi = 32 hu). Here it seems to be appealing to es- 
timate if and where the LOFF phase could appear by 
using the LDA, which should be exact in an infinite sys- 
tem. In order to do this, we calculate at each point R the 
thermodynamic potential Q of a uniform gas with effec- 
tive average chemical potential p e ff(R) = fi — muj 2 R 2 /2, 
assuming a gap of the form A(r) = Aexp(iq • r), and 
minimize £1 with respect to A and q. To be more precise, 
we should have taken a more sophisticated expression for 
the order parameter, but we stress here that we use the 
LDA just to have some indications about what kind of 
behavior one should expect. Let us again look at A(0) 
as a function of the asymmetry (dashed line in Fig. 1). 
We find that LDA gives the BCS solution q = as the 
most favorable solution up to Sfi = 11.9 hw. At that 
asymmetry we find a first-order phase transition (i.e., a 
discontinuity in A(0)) to the LOFF phase with q ~ 17 
which means that the order parameter oscillates with a 
wavelength of ~ 6.2 Iho- This behavior is different from 
the microscopic (BdG) result (solid line in Fig. 1), which 
shows a smooth behavior of A(0). Nevertheless, also in 
the BdG calculation there is a rapid change of A(0) be- 
tween Sfi = 10 kuj (a w 0.25) and Sfi = 11 hw (a w 0.29), 
where the system goes from the BCS-type to the LOFF- 
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FIG. 2: Radial profile of the pairing gap A(r) (in units of fkv) 
for S/j, = 12 huj. The radial coordinate r is expressed in units 
oilho- The microscopic (solid line) and the LDA (dashed line) 
results are plotted. 

type phase, as discussed above. The minimum that one 
observes for the BdG gap at Sfi 10.8 huj and the sub- 
sequent enhancement are due to shell effects which still 
persist even in this large system and which cannot be 
reproduced by the semiclassical LDA calculation. 

In Fig. 2 we plot the radial profile of the pairing field 
A(r) for 5fi = 12 tuo, corresponding to a ~ 0.3, at T = 0. 
The microscopic (solid line) and the LDA (dashed line) 
results are shown. Within LDA, in this case, the LOFF 
phase is more favorable than BCS for all values of r. 
We observe in Fig. 2 that the LDA gap goes abruptly 
to zero at a radius of ~ 2 Z/ lo , which is smaller than the 
LDA wavelength of ~ 6.2 l} lo . Thus, the region where the 
gap is non-zero does not even contain one wavelength of 
the oscillation and therefore the validity of LDA seems to 
be very questionable. As expected from the symmetric 
case [III , the LDA fails to describe the tail of the pairing 
field: The LDA gap goes abruptly to zero while in the 
microscopic case the gap has a smooth profile. We finally 
observe that the microscopic order parameter makes an 
oscillation and that a node is situated at ~ 3 I ho- The 
modulation of the order parameter and the presence of 
a node are signals which indicate that the system is in a 
LOFF-type phase. 

Let us consider now the case of a smaller asymmetry, 
5 fx = 6 hw, corresponding to q « 0.15. We show in 
Fig. 3 the radial profiles of the gap A(r) (top) and of 
the densities (bottom) at T = 0. In the upper panel 
we report the microscopic (BdG) gap (solid line) and 
the LDA result (dashed line). In this case, according to 
the LDA, the BCS phase (q — 0) would be energetically 
preferred in the center of the gas (as we have shown in 
Fig. 1) and up to r — 3.8 Iho, while the LOFF phase 
with q ~ 0.7 would be more favorable in the interval 
3.8 Iho < r < 4.1 Iho- For larger values of r, the LDA 
predicts that the gap should be zero. The wavelength 
of the oscillation of the order parameter in the LOFF 
phase (~ 8.9 Iho) would again be much larger than the 
region where the gap is non zero. Contrary to the LDA, 
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FIG. 3: Radial profiles of the pairing gap A(r) in units of hui 
(top) and of the densities in units of i^j (bottom) for S/j, — 6 
huo as a function of the radial coordinate r in units of Iho- In 
the upper panel the solid and dashed lines correspond to the 
BdG and LDA results, respectively. In the lower panel, the 
solid and the dashed lines refer to the BdG results for the + 
and — densities, respectively. 



the microscopic BdG calculation gives a smooth behavior 
of the order parameter. Near the center, it corresponds 
rather well to the LDA prediction, indicating that the 
pairing is more or less of BCS type. Between r = 4 
Iho and 10 l/ lo , the gap is oscillating (see inset in the 
upper panel of Fig. 3), indicating the appearance of the 
LOFF-type phase. Although within the BdG calculation 
there is no sharp transition from one phase to the other, 
qualitatively it seems that both phases can be present at 
the same time in different regions of the system. 

In the lower panel of Fig. 3 the BdG results for the 
two densities (/?+ and p_) are shown. One observes that 
in the center of the gas the two densities are equal. This 
is coherent with the fact that in the BCS phase at T = 
the LDA always gives p + = p_ if A > 6fi/2, as it is the 
case here. 

All the results shown so far refer to T = 0. However, 
in real experiments with trapped atomic gases the tem- 
perature is always non-zero. Let us therefore raise the 
temperature in the case of asymmetry Sfi = 6 fku in or- 
der to analyze what happens to the gap modulation when 
the temperature is finite. In Fig. 4 we show the order pa- 
rameter in the radial interval from 4 Iho to 10 Iho (where 
we observed an oscillation in the case T = 0) for four 
values of temperature, T = (solid line), T = 0.5 huj/ks 
(dashed line) T = hu/kB (triangles) and T = 2 hiojkB 
(circles). One observes that the oscillation has a smaller 
and smaller amplitude with increasing temperature and 
that it disappears between T = hw /ks and T = 2 Hu/ks- 
Our interpretation of this result is that the critical tem- 
perature of the LOFF-type phase is smaller than the BCS 
critical temperature. Therefore, the LOFF-type phase 
disappears at some temperature between fuv/ks and 2 
hw/ksi while the BCS gap in the central region of the 
gas is still different from zero at T = 2 hw/ks- In this 
case the LDA results (not shown) are very different from 
the BdG ones (the gap is much too large), as one could ex- 
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FIG. 4: Radial profile of the pairing gap A(r) (in units of fkv) 
for <5;U = 6 hui. The radial coordinate r is expressed in units 
of Iho- Results at T — (solid line), T = 0.5 huo/ks (dashed 
line) T = huj/ks (triangles) and T = 2 huj/kB (circles) are 
reported. 

pect from the fact that already in the symmetric case the 
agreement between LDA and BdG calculations becomes 
quite poor at finite temperature 0. Nevertheless, qual- 
itatively the LDA gives again the right indication: Also 
within LDA the LOFF phase is absent at T — 2 hw/ks- 
Instead, at that temperature the Sarma phase becomes 
more favorable in certain regions within the LDA: This 
phase is characterized by zero momentum Cooper pairs, 
a gap A smaller than Sfi/2, different densities p+ and p- 
and typical occupation number distributions as shown in 
Ref. 0. 

To summarize, we have solved the BdG equations for 
an atomic Fermi gas with different populations of two hy- 



perfine states. It is well-known that an increasing asym- 
metry of the populations renders BCS pairing difficult, 
and non-standard pairing mechanisms become possible. 
In this article we showed that the BdG formalism au- 
tomatically includes such non-standard pairing mecha- 
nisms through the non-diagonal matrix elements of the 
gap. For example, in the case of a small system, we found 
that the usual pairing disappears at a certain asymme- 
try, but when the asymmetry is strong enough such that 
the single-particle energies of states with opposite spin 
and different principal quantum numbers start to match, 
pairing becomes again possible, but now with an oscillat- 
ing order parameter. This is very similar to the LOFF 
phase introduced for the case of a homogeneous system. 
In the case of a larger system, there is no longer a sharp 
separation between the BCS pairing and the LOFF-type 
pairing: As a function of asymmetry, but also as a func- 
tion of the distance from the center of the trap, the sys- 
tem undergoes smooth transitions from one kind of pair- 
ing to the other. This result is qualitatively different from 
that obtained with LDA calculations, where the transi- 
tion between the BCS and the LOFF phase is a first 
order phase transition. We also observe that even a sys- 
tem containing 17000 atoms is still much too small for 
the LDA to be applicable, since the wavelength of the 
LOFF oscillations is of the same order of magnitude as 
the whole system. Finally we looked at the temperature 
dependence of the LOFF-type phase. We observe that it 
disappears already at temperatures where the BCS phase 
is still present. This, of course, can be a problem if one 
tries to observe the LOFF phase in experiments. 

We acknowledge discussions with M. Baldo, F. Catal- 
iotti and A. Sedrakian. 
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